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We report density functional calculations of the electronic structure, Fermi surface, phonon spec- 
trum, magnetism and electron-phonon coupling for the superconducting phase FeSe, as well as the 
related compounds FeS and FeTe. We find that the Fermi surface structure of these compounds is 
very similar to that of the Fe-As based superconductors, with cylindrical electron sections at the 
zone corner, cylindrical hole surface sections, and depending on the compound, other small hole 
sections at the zone center. As in the Fe-As based materials, these surfaces are separated by a 2D 
nesting vector at (7r,7r). The density of states, nesting and Fermi surface size increase going from 
FeSe to FeTe. Both FeSe and FeTe show spin density wave ground states, while FeS is close to 
an instability. In a scenario where superconductivity is mediated by spin fluctuations at the SDW 
nesting vector, the strongest superconductor in this series would be doped FeTe. 

PACS numbers: 74.25.Jb,74.25.Kc,74.70.Dd,71.18.-|-y 



I. INTRODUCTION 

Superconductivity was recently reported in a-FeSei-^^, 
with critical temperature Tc ^ SKi This is of inter- 
est both because of the fact that Fe containing su- 
perconductors are unusual and because this material 
shares in common square planar sheets of tctrahcdrally 
coordinated Fc with the Fc-As based high tempera- 
ture superconductorsi^'^'^'^'i a-FeSe occurs in the PbO 
structure. This consists, as mentioned, of Fe square pla- 
nar sheets, with Se atoms forming distorted tetrahedra 
around the Fe very similar to the structure of the FeAs 
planes in LaFeAsO, BaFe2As2 and LiFeAs, which are pro- 
totypes of the known families of Fe-As based high-Tc 
superconductors i^ii^ 

Since the reported Tc = 8K of doped a-FeSe is mod- 
est, it is important first of all to establish the relationship 
between this material and the Fe-As based superconduc- 
tors. We note that LaNiPO is also a superconductor and 
shares the crystal structure of LaFeAsO r^^ but that it 
is apparently quite different electronically and can be un- 
derstood in terms of standard electron-phonon theory^i^ 
unlike the Fe-As based phasesJ^ii^ 

Here we report density functional calculations that 
show a- FeSe and the other known Fe based chalcogenides 
in this structure to be very similar to that of the Fc-As 
based superconductors li^ In particular the Fermi surface 
consists of small heavy hole cylinders near the zone cen- 
ter and lighter compensating electron cylinders around 
the zone corner. We show that the stoichiometric com- 
pounds are either very close to a spin density wave in- 
stability (FeS) or have an itinerant spin density wave in- 
stability without doping (FeSe and FeTe) similar to the 
Fe-As superconductors iii We predict that this itinerant 
nesting driven magnetic state is strongest in FcTc and in 



addition that FeTe has the largest Fermi surface of the 
three compounds. Calculations of the electron-phonon 
coupling show that doped FeSe is not an electron-phonon 
superconductor, similar to what was found for the Fe-As 
phases. Within a spin-fluctuation driven picture of super- 
conductivity the results indicate that FeTe with doping 
is a likely higher temperature superconductor. 



II. FIRST PRINCIPLES METHODS AND 
STRUCTURE 



Our calculations of the electronic structure and mag- 
netic properties were performed within the local den- 
sity approximation with the general potential linearized 
augmented planewavc (LAPW) methodfi^ including lo- 
cal orbitalsfi^ similar to our previous calculations for 
the Fe-As based superconductors J^ii^i2£ We used LAPW 
spheres of radius 2.1 ao for Fe, Se and Te and 1.9 oq for 
S. These compounds occur in a simple tetragonal struc- 
ture with one internal parameter, zx corresponding to 
the chalcogen height above the Fe square plane. The ex- 
perimental lattice parametersii2ii^ were employed and 
we relaxed the chalcogen height via energy minimiza- 
tion. The structural parameters used and some results 
arc presented in Table IH The electron-phonon coupling 
and phonon dispersions were on the other hand obtained 
using linear response, again with the experimental lat- 
tice parameters, with the Quantum Espresso code^i^ 
within the generalized gradient approximation of Perdew, 
Burke and Ernzcrho^ as described for LaFeAsO and 
LaNiPOJ^ii^ 



TABLE I: Structural parameters and magnetic properties of 
PbO-structure FeX. The lattice parameters are from exper- 
imental data, while the internal chalcogen structural param- 
eter, zx is from LDA structure minimization, msow is the 
spin moment within the Fe LAPW sphere (radius 2.1 ao) for 
the SDW state and Esdw is the energy per Fe of this state 
relative to the non-spin-polarized state in meV/Fe. N{Ep) 
is the density of states at the Fermi energy in the non-spin- 
polarized band structure in eV~^ on a per Fe both spins basis. 

a(A) c(A) z^ N{Ef) msDw{fJ.B) Esdw 



FeS 3.6735 
FeSe 3.765 
FeTe 3.8215 



5.0328 
5.518 
6.2695 



0.2243 
0.2343 
0.2496 



1.35 
0.95 
1.83 



0.00 
0.65 
1.28 
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III. ELECTRONIC STRUCTURE, PHONONS 
AND MAGNETISM 

Our main results for the electronic structure are given 
in Figs. [U [2] and [31 which show the non-spin-polarized 
band structures, electronic densities of states and Fermi 
surfaces of FeS, FeSe and FeTe. The calculated values 
of N{Ep) are given in Table HI The phonon dispersions 
of FeSe are shown in Fig. [H One interesting feature of 
the phonon dispersion is that they have little dispersion 
in the kz and e.g. are very flat along the tetragonal T-Z 
direction. This presumably reflects anion-anion repul- 
sion, which leads to long bonds between the FeSe layers, 
the result is that there may be an easy cleavage plane 
between the Se ions, which may facilitate preparation 
of clean surfaces for experiments such as photoelectron 
spectroscopy. 

The phonon density of states, G{uj) and electron- 
phonon spectral function, a'^F{LLi) are given in Fig. \5\ 
The electron-phonon coupling constant for FeSe as ob- 
tained in linear response is A = 0.17 with a;;og=113 cm~^. 
No superconductivity at any temperature even approach- 
ing IK results with these values within standard electron- 
phonon theory even if very low values of the Coulomb 
parameter, e.g. ^* = 0.10 are used in the Allen-Dynes 
equation. This is similar to what was found previously 
for LaFeAsO ji^ii^ Therefore we conclude that FeSe is not 
a conventional electron-phonon superconductor. 

Turning to the electronic structure we find a strong 
qualitative similarity between these materials and the 
FeAs-based superconductors. In particular, we find these 
to be low carrier density metals, with high density of 
states. This arises from band structures that are closely 
related to those of the FeAs materials. The chalcogen p 
states lie well below the Fermi level and are only modestly 
hybridized with the Fe d states as may be seen from the 
projected DOS (Fig. [2]). Thus the electronic structure 
near the Fermi energy derives from metallic Fe^"*" lay- 
ers, with direct Fe-Fe interactions. These are embedded 
inside a largely ionic background which imposes a com- 
peting tetrahedral crystal field. As in the Fe-As based 
materials, there is a pseudogap at an electron count of 6 
d electrons per Fe, and Ep lies near the bottom of this 
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FIG. 1; Band structures of FeS (top), FeSe (middle) and FeTe 
(bottom) from non-spin-polarized calculations with the LDA 
relaxed X heights. 



pseudogap. We emphasize that this is not the position 
of a tetrahedral crystal field gap, which would be at 4 
electrons, and emphasizes the fact that Fe chalcogen hy- 
bridization is not strong compared with the Fe-Fe interac- 
tions. This explains the similarity of the electronic struc- 
ture to that of the FeAs-based materials, which were also 
found to be substantially ionic in similar calculations^^ 
These band structures yield two intersecting elliptical 
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FIG. 3: (color online) LDA Fermi surface of FeS, FeSe and 
FeTe from non-spin-polarized calculations with the LDA re- 
laxed X heights. The corners are F points. 




FIG. 2: (color online) Electronic DOS and projection onto 
the LAPW Fe and chalcogen spheres indicating the Fe d and 
chalcogen p contributions for FeS (top), FeSe (middle) and 
FeTe (bottom) as in Fig. [T] 
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FIG. 4: Calculated GGA phonon dispersions of non-spin- 
polarized FeSe. 



cylindrical electron Fermi surfaces at the zone corner in 
all three materials. These are compensated by lower ve- 
locity hole sections at the zone center - two concentric 
hole cylinders, and in the case of FeS and FeTe a small 
closed hole section inside the inner cylinder. This is qual- 
itatively very similar to the FeAs-based materials. It is 
important to note that cylinders at the zone center and 



zone corner, if they are the same size, would yield strong 
nesting peaked at the 2D (7r,7r) point. This will lead 
in general to enhanced spin fluctuations at the nesting 
vector, and if sufficiently strong will cause a spin den- 
sity wave. In fact, an ordered SDW was found both 
in first principles calculations and in experimental stud- 
ies for LaFeAsO and many other undoped FeAs-based 
compounds] ^^'^^i^^'^^i^^'^^i'^°''^^ There is a clear competi- 
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FIG. 5: Calculated GGA phonon density of state G{uj) and 
electron-phonon spectral function a^F{u>) for FeSe. 



tion between the SDW and the superconducting state 
in that superconducting samples generally do not show 
the SDW, while samples with the SDW transition, gener- 
ally do not show superconductivity. The ground state is 
an antiferroniagnetic cell doubled along the [11] in plane 
direction to yield linear chains of nearest neighbor like 
spin Fe atoms arranged antiferromagnetically, although 
the values of the moments are dependent on details, es- 
pecially the As height above the Fc plane i^2iSiMi2^i2&i21 



IV. DISCUSSION 

Turning to the trends, the size of the pseudogap is ap- 
proximately the same in FeSe and FeS, but is significantly 
smaller in FeTe. Specifically, there is a greater overlap 
between the hole and electron bands in the latter com- 
pound. This leads to larger Fermi surfaces. The value of 
N{Ep)=:l.83 eV~^ is also highest in this compound, al- 
though it is still lower than the 2.62 eV~^ that is obtained 
for LaFeAsO by the same approach. We used a super- 
cell approach to investigate the SDW with the chalco- 
gen heights fixed to the values calculated by non-spin- 
polarized energy minimization. We find instabilities for 
FeSe and FeTe, but not for FeS. The spin moments and 
the energy of the SDW relative to the non-spin-polarized 
state are given in Table [H As may be seen the SDW is 
considerably stronger in FeTe than in FeSe, with an en- 
ergy gain of 47 meV/Fe and a spin- moment of 1.3 fiB- 
The corresponding values for LaFcAsO, calculated in the 
same way are Esdw=^^ meV/Fe and msow = 1-0 Hb 
in a 2.1 uq Fe sphere. Besides the SDW we find a bor- 
derline ferromagnetic tendency in FeTe, without doping 
when the SDW is not allowed, similar to LaFeAsO.— 
We do not find ferromagnetic instabilities in either FeSe 
or FeS, consistent with the lower values of N{Ep) in 
those materials. The sensitivity of the moment size to 



the ordering pattern underscores that fact that these are 
itinerant magnetic systems in the LDA, as opposed to 
local moment magnets. This means that magnetic or- 
dering is driven by electrons at and near the Fermi sur- 
face. On the other hand this is not to say that spin fluc- 
tuations are weak in the paramagnetic state above the 
SDW ordering temperature or the paramagnetic state as 
realized by doping. In fact, as noted, several authors 
have found that the As height in the Fe-As compounds 
is strongly coupled to magnetism and so strong spin fluc- 
tuations in the paramagnetic state would help rational- 
ize the underestimated As height in non-spin-polarized 
LDA calculations42iSi2ii^i22 In fact, there is evidence 
for strong spin fluctuations in the normal state of the 
FcAs compounds, e.g. from temperature dependent re- 
sistivity data indicating strong scattering above the SDW 
ordering temperature (note a drop in resistivity below the 
SDW ordering temperature even though the carrier den- 
sity is strongly reduced) ^^ as well as from spectroscopy i^ 
The size of the effects observed implies that these fluctu- 
ations should have large amplitudes and therefore unlike 
the SDW should be rather diffuse in q-space, which might 
make them hard to directly observers Furthermore, the 
strong transport signatures especially the enhanced resis- 
tivity above the ordering temperature imply substantial 
coupling between spin flucuations and electrons at the 
Fermi energy, which is indicative of the itinerant nature 
of the magnetism. 

We emphasize also that our results in Tablc|T]are at the 
LDA relaxed atomic positions for the non-spin-polarized 
states. These systems become more magnetic as the 
chalcogen height is raised. In the case of FeS a recent 
refinement is available, and gives 2:3=0. 2602i^ This puts 
the S ions 0.18 A higher than in the LDA structure. 
With this value we find a stable SDW for FeS, with a 
moment of 1.2 /Lts/Fe. This is the same trend as in the 
Fe-As based superconducting materials. The result may 
be taken as an indication that in fact FeS may have a 
spin density wave as well and at least that there will be 
strong spin fluctuations in FeS as well as the other PbO- 
structure Fc chalcogenides. 

As mentioned, cylindrical Fermi surface sections of 
equal volume will be nested with nesting vector equal 
to the separation of the centers of these cylinders. This 
nesting can be reduced by imperfect matches in shape, 
three dimensionality, and size mismatch. Size mismatch 
can arise both from additional Fermi surface sections, 
such as the extra small hole sections obtained in FeS and 
FeTe, but not FeSe, or from doping. In particular, elec- 
tron doping will reduce the size of the hole sections and 
increase the size of the electron sections consistent with 
Luttinger's theorem. 

The mechanism for superconductivity in the Fe-As 
based phases has not yet been established. Nonethe- 
less, there are indications that magnetism is associated 
with superconductivity. These include the modest elec- 
tron phonon couplings in the materials, the proximity to 
magnetism and the phase diagrams which show an as- 



sociation between the SDW and superconductivity. We 
discuss our results within a general spin fluctuation me- 
diated framework fi^^i^i^ In general an itinerant SDW 
instability arises from a divergence of the real part of the 
susceptibility x(ct) ^t a specific wavevcctor q. Supercon- 
ducting pairing is also associated with the real part of 
x(q) through an integral over the Fermi surface. x(q) 
for q connecting different parts of the Fermi surface can 
contribute to pairing. As a result, when the Fermi sur- 
faces are small and disconnected as in these Fe based 
materials, the fluctuations associated with nesting will 
provide substantial interband pairing between the elec- 
tron and hole sections, but will not provide substantial 
intraband pairingp^^^ As the nesting is reduced, e.g. by 
doping, x(q) will spread out while the peak value will 
be reduced; consistent with the destruction of the SDW. 
For circular cylinders of radii differing by Sq, x(q) '^iU 
show a plateau of high x(ct) around the nesting vector 
with width 25q (this will persist until the radii differ by 
a factor of two at which point the center of the plataeu 
will dip). This means that the two cylinders will still be 
connected by spin fluctuations associated with their now 
reduced nesting, and therefore even though the SDW will 
be suppressed by the reduction in the maximum value of 
X the associated spin fluctuations can still provide super- 
conducting pairing. We emphasize that in this general 
framework the same parts of the Fermi surface are af- 
fected by the SDW and by the superconducting pairing. 
Thus these two Fermi surface instabilities compete for 
the same Fermi surface and therefore that there should 
be at best little coexistence of these two orders. 



gen atoms in these compounds. As in the arsenides, we 
find that the electron-phonon coupling cannot explain 
the superconductivity, and furthermore that these com- 
pounds display itinerant magnetism. These results imply 
a similar superconducting nature for the Fe-As phases 
and FeSe. 

The trend that we find in going from FeSe to FcTe 
is interesting in this context. In particular we find quite 
cylindrical Fermi surfaces and an SDW instability in both 
compounds. However, the strength of the SDW is sub- 
stantially higher in FeTc as is the size of the Fermi sur- 
face. Within the general framework discussed above, 
FeTe would be expected to have stronger pairing, and 
therefore higher Tc than FeSe assuming that the same 
mechanism applies in both materials, that both mate- 
rials can be chemically doped to the optimum carrier 
density and that competing instabilities do not prevent 
superconductivity in that case. It will be of interest to 
experimentally probe the similarities of FcSe with those 
of the Fc-As phases and to search for superconductivity 
in doped FeTe and in the alloy Fe(Sc,Te). 
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SUMMARY AND CONCLUSIONS 



VI. NOTE ADDED 



We report electronic structures, magnetic properties 
and electron-phonon calculations for FeX, X=S,Se,Te. 
We find strong similarities to the Fe-As based supercon- 
ductors, reflecting the ionic nature of the As and chalco- 



Mizuguchi and co-workers^ recently reported observa- 
tion of superconductivity at 27K in FeSe under pressure. 
This is consistent with the conclusion here regarding the 
relationship between FeSe and the Fe-As based materials. 
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